Microscopic theories of neutrino-^^C reactions 



C. Volpe'^, N.Auerbach^ G.Colo^ T.Suzuki'', N. Van Giai" 
Groupe de Physique Theorique, Institut de Physique Nucleaire, F-91406 Orsay Cedex, France 
School of Physics and Astronomy, Tel Aviv University, Tel Aviv 69978, Israel 
'^^ Dipartimento di Fisica, Universita degli Studi, via Celoria 16, 20133 Milano, Italia 
Departement of Physics, College of Humanities and Science, Nihon University, Sakurajosui 3-25-40, Setagaya-ku, Tokyo 

156, Japan 

In view of the recent experiments on neutrino oscillations performed by the LSND and KARMEN 
collaborations as well as of future experiments, we present new theoretical results of the flux averaged 
^'^C{ue,e-y'^N and '^C(i/^, m" )^^A cross sections. The approaches used are charge-exchange RPA, 
charge-exchange RPA among quasi-particles (QRPA) and the Shell Model. With a large-scale shell 
model calculation the exclusive cross sections are in nice agreement with the experimental values 
for both reactions. The inclusive cross section for V/j, coming from the decay-in-flight of is 
15.2 X 10~*° cm? (when Hartree-Fock wavefunctions are used), to be compared to the experimental 
value of 12.4 ± 0.3 ± 1.8 x 10"*" cm^, while the one due to i^e coming from the decay-at-rest of fi'^ is 
16.4 X 10^*^ cm^ which agrees within experimental error bars with the measured values. The shell 
model prediction for the decay-in-flight neutrino cross section is reduced compared to the RPA one, 
namely 19.2 x 10"''° cm^. This is mainly due to the different kind of correlations taken into account 
in the calculation of the spin modes (in particular, because of the quenching in the 1^ channel) and 
partially due to the shell-model configuration basis which is not large enough, as we show using 
arguments based on sum-rules. Results for exclusive and inclusive muon capture rates and beta 
decay are given and are close to the experimental findings. 



I. INTRODUCTION 



For many years, weak processes in nuclei such as beta decay and muon capture have been studied on the one hand to 
deepen our knowledge of the weak interaction in nuclei and on the other hand to yield information on nuclear structure. 
Both aspects are important when one considers other weak processes in nuclei, namely reactions induced by scattered 
neutrinos. A description of the latter processes is not only important in our attempts to better understand the nature 
of such reactions but also it has significant practical importance in current experimental studies of neutrinos. In 
fact, nuclei are often used as neutrino detectors so that the knowledge of reactions induced by neutrinos on nuclei 
becomes a crucial step for the interpretation of experiments on neutrinos, such as the ones aiming to go beyond the 
Standard Model looking for neutrino oscillations and masses or those measuring solar neutrinos to test the Standard 
Solar Model. 

A clear example is given by the recent experiments performed both by the LSND and the KARMEN collaborations, 
looking for i/^ —^ i^^i Vf. [||J^] or t'^ ^ @ oscillations with neutrinos produced by accelerators. The 

detectors used in these measurements are mainly composed of protons and ^^C. Reactions of neutrinos on this nucleus 
are used to check neutrino fluxes and efficiencies In the v^, — > disappearance experiment of Ref. |^ with v^. 
coming from the Decay- At-Rest (DAR) of the charged-current (CC) reaction v^^ -|-^^ C ^ e~ -l-'^^ ^g.s. is used to 
detect neutrinos. In Ref. the appearance experiment with v^^ coming from the Decay-In-Flight (DIP) of 

Tr"*", ^^C is used to detect neutrinos via the CC reaction Vf. -|-^^ C -f N . In both cases, the extracted oscillation 

probabilities rely directly on the knowledge of the cross section for these reactions. The exclusive cross section for 
j/g C — > e~ -f Ng.s. (where ^"^N is left in the ground state) has been measured by different collaborations 
and its value can be obtained in a model independent way by using form factors deduced from the measurements of 
related weak processes such as beta decay and muon capture as it is done in the Elementary Particle Theory (EPT) 
pO| . Concerning the inclusive cross section for Ve -f^^ C ^ e~ -|-^^ TV (where ^^iV is left either in the ground state 
or in an excited state), because of the universality of the weak interaction we expect this reaction to be described 
by the same effective hamiltonian as the reaction C ^ TV with i^^ coming from the DIP of 7r+. A 

problem concerning this reaction has emerged which has been extensively investigated recently, nainely the theoretical 
cross section overestimates the experimental value by about 50% within charge-exchange RPA |0,|2| or by about 
30 — 40% in An attempt has also been made within an extension of the EPT lO], but it is based on several 

assumptions which have not been tested yet jl^ . In a recent shell model calculation [|l6| the value of the cross section 
is reduced in apparent agreement with the experimental value. The reaction cross section with neutrinos coming from 
the DIP of 7r+ has been measured only by the LSND collaboration up to now [6|JT7|] . Because the same detector used 
for these measurements has been used to measure the neutrino oscillations, it is important to know what we can say 



1 



about these reaction cross section from the theoretical point of view, keeping in mind that every nuclear structure 
model necessarily contains approximations. 

The weak interaction in nuclei being well known, an accurate prediction of these cross sections is a challenging 
problem from the nuclear structure point of view. In fact, these observables rely on transition densities which can 
be obtained only in calculations which take into account the different aspects of the structure of the nuclei involved 
in these reactions. First, it is known for a very long time jl^ that ^^C is not a true closed sub-shell nucleus or, 
in other words, that the ground state wave function can be only described by an intermediate coupling scheme and 
contains configuration mixing, including deformed components. Second, if on one hand in the DAR experiments the 
neutrinos have impinging energies of the order of several tens of MeV, in the DIF experiment the energy goes up 
to around 300 MeV. As a consequence, the DAR reaction cross section is dominated by the Gamow- Teller (GT) 
transition to the ground state of ^^iV for which information on the transition probabilities can be obtained from other 
related weak processes like ^'^N{f3+)'^'^C, ^'^B{/3-y^C and ^'^C{iJ,~y^B. On the contrary, in the DIF reaction the 
energy and momentum transferred to the nucleus is quite large so that ^'^N can be left in an excited state of several 
tens of MeV, that is, in and above the giant resonance region. If the transition to ^'^Ngs represents 2/3 of the total 
reaction cross section in the DAR case, the cross section given by transitions to the excited states is 200 times larger 
than the transition to the ^^iVg^, if the neutrinos come from the DIF of pions. This means that to have an accurate 
prediction for these cross sections we need nuclear structure models which are not only capable of taking into account 
configuration mixing in the ground state wave function of ^^C but also of describing high-lying states in ^'^N. This is 
a quite difficult task within the present nuclear structure models. 

The microscopic theoretical approaches used so far are either the charge-exchange Random-Phase- Approximation 
(we will call it from now on simply RPA) or the shell model . The former includes only partially configuration 
mixing in the ground state while it can easily include high-lying one particle-one hole (Ip-lh) configurations. On the 
other hand the shell model can give a good description of the ground state wave function whereas the prediction of 
high-lying states requires a large model-space which may be very difficult to treat numerically. 

In this paper we try to improve the existing calculations within the two microscopic approaches just mentioned. 
First, we use a charge- exchange RPA approach applied to quasi-particles (we will call it from now on simply QRPA), in 
order to improve the poor description of the ground state of ^^C within RPA. The configuration mixing is introduced 
by including "ad /loc" pairing correlations in this nucleus. This is done to see if the configuration mixing that is 
missing in the RPA ground state wave function can be at the origin of the discrepancy between the measured reaction 
cross section and the theoretical predictions as it was first suggested in It was shown in that the inclusion 
of fractional occupancies within the RPA approach reduces the exclusive cross section from a factor 3-4 to 50% 
discrepancy. In Ref. |2|] an RPA calculation with partial occupancies was performed and it led to similar results for 
the exclusive cross section while the inclusive cross section was shown to decrease by a few percent. However, the 
authors of flj,^ did not show the results of a full QRPA calculation as we do here. 

Next, we perform shell model calculations within a large model space. Our space is larger than that actually used 
in jl^ where the results in the same space we use are obtained by extrapolation only. We will compare our results to 
the experimental findings and to the other microscopic theoretical predictions We will conclude 

by summarizing the present status of the problem. 

In section II we briefiy review the general theory describing neutrino scattering on nuclei, muon capture and beta 
decay. In section III we present the essential features of the microscopic models used in our calculations, i.e. the 
QRPA and the shell model. In section IV results from the above mentioned weak processes are presented focusing in 
particular on the cross sections for the reactions vi C I N {I — e, /i), both inclusive and exclusive, with z/g 
coming from DAR of /x"*" and with coming from the DIF of 7r+. Conclusions are drawn in section V. 



II. GENERAL THEORY 



The theoretical framework to study nuclear responses to weak probes is discussed extensively in the literature. A 
detailed description can be found in | |2^ . Here we will just present the main ingredients of the calculation necessary 
for the discussion. 

The general expression for the cross section of the reaction vi C I N (l = e, ii) is 




where Ef (Ei) is the energy of the final (initial) nuclear state, E^, (p^) is the incident neutrino energy (momentum) 
and El (pi) is the outgoing lepton energy (momentum). The effective single-particle hamiltonian Heg is derived by 
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carrying out the Foldy-Wouthuysen (FW) transformation and retaining terms up to 0{{\q\/M)^) (M is the nucleon 
mass and q is the momentum transfer), since the momentum transfer involved in nucfear scattering of neutrinos 
produced by acceferators can be large. The expression for Hf,ff can be found in [^. If the nuclear recoil effects are 
ignored, we have 



a = —cos^Oc'^PiEi j d{cosd)Mf3, 



(2) 



where G cos 9c is the weak coupling constant, 9 is the angle between the directions of the incident neutrino and the 
outgoing lepton and Mp is given by 

Mp = MF\{f\l\i)\^ + MGo^|(/|a|i)|2 + Mg2A (3) 
where the squared nuclear matrix elements are 

A-TT 

\{f\m' - ^2J, + 1) ^ l^-^^'" ^ t+{k)Mqrk)Y,{r,)U)\\ (4) 

A'jr 

= T^T-T-) T.\W\ E t+{k)ji{qn)[Yi{f,) x J,)p, (5) 



11'2\(112] 4n 



A = (-)^ y (-i)'/2-''/2+/f72mv2FTT| ,,,, 

° IV K V T V -r y^QQQj (2J, + 1) 

X (j^||^i+(fcU(grfe)[rKr"fc)xcT](^)||j.) 

k 

X {Jf\\J2t+{k%,{qrk-)[Yi,if,,) X a](^)||J,)t. (6) 

k' 

where k labels the space and spin-isospin coordinates of the fc-th nucleon, I, I' are the orbital angular momenta and K 
is the total angular momentum of the transition operators. The coefficients Mp, Mqq and Mq2 appearing in (H) 
depend on the momentum transferred to the nucleus [q = (q, iqa) = pi — Pv) and the standard nucleon form factors 
/v(?^)7/4(9^),/w(9^),/p(9^)- Second-class current form factors are ignored. 

A correction to (|^) must be introduced to account for the distortion of the outgoing lepton wave function due to 
the Coulomb field of the daughter nucleus. For reactions on ^'^C with neutrinos from the DAR of /i^, the quantity 
Pi Ra {Ra is the radius of the nucleus) is of the order of 0.5. In this case, the cross section (||) may be multiplied by 
the Fermi function F{Zf, Ei) [Q, where Zf is the charge of the daughter nucleus and Ei the energy of the charged 
lepton. When the neutrinos come from the DIF of 7r+, the outgoing muons havepii?^ > 0.5. With relativistic leptons, 
the effect due to the Coulomb field may be included by using the "Effective Momentum Approximation" (EMA) . 
In this approximation the lepton energy and momentum are modified by a constant electrostatic potential within the 
nucleus, i.e. E^eff ^ Ei ~ V{Q) and pi^ejf = (-^L// " m'^f/'^ with ^(0) = -2>Zfa/2R. The cross section @ has 
then to be multiplied by a factor pi^ef fEi^ef f / (piEi). 

To obtain the flux-averaged cross sections (a) f that can be compared with the experimental data, the energy 
dependent cross section (^) has to be folded with the normalized neutrino flux f{E^) (depending on the neutrino 
source used) 



{a)f = J dEMEu)f{E,), (7) 

where 

^^""'^ = lEjKfiEir 

f{E^) being the initial flux and Eq the threshold energy. 
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Closely related to this neutrino capture reaction is the capture of a negative muon bound in an atomic orbit, 
/i" + (A, Z) (A, Z — 1)* -\- Uf^. In the IS-capture the inclusive rate Ac is given by 

9 

771 

Ac = V^l^isHG'yM^ + GlMl + {Gl - 2GpGa)MI]. (9) 

if we neglect the recoil term which represents a correction of a few percent . The function (pis is the muon IS-bound 
state wave function evaluated at the origin, i.e., |0isp = R{Zam'_f/n, R being a reduction factor accounting for the 
finite size of the nuclear charge distribution {R = 0.86 for ^^C) and m' the muon reduced mass. The constants 
Gv,Ga,Gp ||2^ are the "effective coupling constants" which depend only slightly on the neutrino momentum p^. 
This can be simply obtained from the energy and momentum conservation, = m' — (to„ — nip) — — i?/;, 
where mn,mp are the neutron and proton masses, is the binding energy of the muon in the IS* orbit and Efi 

is the nuclear excitation energy measured with respect to the parent nucleus ground state. The capture rate can be 
factorized as in (||) if we neglect the dependence of the coupling constants on p,^. The square of the vector, axial- vector 
and pseudoscalar matrix elements are 



I f k 



m 

,K f 



=47r ^(2/^ + l)^(^)^|(Jyl5]t+(A:)0p(i^)|J,)P, (12) 

l,l',K f ^ k 



where 



Op{K) = ^2K^^f iP'^^^)iYiih) X + ^^&-i)^2j,'(p.r,)[y,(f,) x . (13) 

with I ^ K — 1 and I' ^ K + 1. The /3-decay corresponds to the limit of zero momentum transfer of the transition 
probabilities in (||), that is 

M° ^ fm\{f\m' + /i(o)^i(/kK>p (14) 



where 



I(/|1N>P = ^;^I(^/IIE^+WII-^^)I' (15) 



\{fW\i)\' = ^J^liJfW E t+{kMk)U}\ (16) 

k 

and the transition operators are of the usual Fermi or Gamow- Teller type. The ft value is given by 

^ [G^cos-^Ocml) Mj' ^ ' 

The neutrino reaction cross section (^) , the muon capture rate (||) and the ft value for the beta decay (Q) depend 
on the wavefunctions of the initial and final nuclear states involved in these processes. The microscopic models used 
in this work to evaluate these wavefunctions and the corresponding transition probabilities are described in the next 
section. 
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III. MICROSCOPIC MODELS 



A. RPA and QRPA 

Transition matrix elements of the type entering in eqs. flplJl^) can be calculated within the framework of RPA 
or QRPA. In the present work, the starting point is a Hartree-Fock (HF) calculation of the ground-state of '^^C, 
performed in coordinate space by using the Skyrme-type effective interactions SIII ||2^ and SGII (29|. The SGII force 
was built with the purpose of obtaining a proper description of spin-isospin nuclear properties. The HF solution 
determines the mean- field and single-particle (s.p.) occupied levels of ^^C which has a closed-subshell structure in this 
description. The unoccupied levels of ^^C are obtained by diagonalizing the HF mean-field using a harmonic oscillator 
basis. Therefore, the continuum part of the s.p. spectrum is discretized and discrete particle-hole (ph) configurations 
coupled to J"^ are used as a basis in order to cast the RPA equations in the matrix form. The details of this procedure 
can be found in . This RPA calculation is self-consistent since the residual interaction among ph states is derived 
from the same Skyrme force used to produce the mean field. 

To go beyond the closed-subshell approximation for the ^^C ground state, pairing correlations are taken into account 
in the HF-I-BCS approximation. Constant pairing gaps Ap and A„ for protons and neutrons are introduced and are 
set at 4.5 MeV. A large pairing gap is unrealistic for states far from the Fermi surface, and an energy cut-off is 
required, such that the states above this cut-off have A = 0. The cut-off is set at the 2si/2 state. On top of the 
HF-I-BCS calculation, the QRPA matrix equations can be written with a procedure which parallels what was described 
above, with the two-quasiparticle (2qp) configurations replacing the ph ones. We do not present here the details of the 
QRPA formalism, which is found in the literature (see, e.g., Rcf. [Q). We simply note that the particle-particle matrix 
elements are here renormalized by means of a parameter that has been chosen to be smaller than 1 (typically 0.7) 
to avoid the well-known ground-state instabilities. 

For the multipolarities studied, it is found that the space used for the RPA (QRPA) calculations satisfies well (up 
to a few percent) the energy-weighted sum rule (EWSR) associated with the operators of the type 

Y,t+{k)riY,{v^), ^t+(fc)ri[l^(ffe)xa](-') (18) 

k k 

which are the small-g limit of those defined in the preceding section. 

For each multipolarity, every eigenstate of the RPA or QRPA equations is characterized by its and amplitudes 
and the transition matrix element for a generic operator 0{k) is written as 

{Jf\Y,d{k)\.h)=Y.^a\Y,d{k)mxipU^vp + YlpV^up), (19) 

k a,f3 k 

where a and P label a given ph or 2qp states, u and v are the BCS occupation amplitudes (which reduce to 1 and 
in the HF-RPA case) and (a| J^k single-particle matrix elements. With A — 4.5 MeV, the occupation 

probabilities for the lp3/2 and lpi/2 states are Vp^^^ = 0.84 and Wp^^^ ~ ^-^^ respectively. The latter one is smaller 
than the one used in [|l2|;^ . 



B. Shell Model 



We also evaluate transition matrix elements for neutrino and muon capture reactions using a large shell model space. 
It is desirable to use extended spaces as much as possible in order to treat both inclusive and exclusive reactions. 
Here, we take the Os-Op-lsOd-lpOf shell model space and include configurations up to "ihuo excitations for negative 
parity states and up to 2fujj excitations for positive parity states. No *He core is assumed in the present calculations. 
The configurations taken for positive parity states in A=12 nuclei are 

+ (0s)3(0p)«(ls0d)i + (0s)2(0p)i*^ (20) 

and those for negative parity states are 

(Os)'*(Op)^(lsOd)i + {Qsf{Qpf + {fdsf{Qpf{ls{)df{lpQfY 
+ [Osf{Opf{lsMf + {Osf{Opf{lpQfY 

+ [Qsf{Qp)\lsMf + [Qsf{Qpf{lsQdY (21) 
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In the shell model calculations of ref. [|T6[ of the present problem, computations were carried out with configurations 
up to Ihuj excitation for negative parity states and up to 2hu} excitation for positive parity states. 

The spurious center-of-mass states are eliminated here by using the method of Lawson . The spurious states 
are pushed up to higher energies with the addition of a fictitious term in the hamiltonian which acts only on the 
center-of-mass excitations. The number of spurious states removed is ^ 1000 — 4000 for each negative parity multipole 
(0~, . . . , 4^) and ~ 100 — 400 for each positive parity multipole (0+, . . . , 6+). The number of states (without the 
spurious center-of-mass components) amounts to ~ 2800— 10600 for the negative parity multipoles and ^ 250— 1500 
for the positive parity multipoles. The number of states is very large for the negative parity multipoles compared to 



the work in ref. |16|, in which the number was ~ 50 — 150 and included up to Ihuj excitation only. 

We adopt here the effective interaction of Warburton and Brown for use in the present Os- Op- IsOd-lpOf model 
space, and we use the set WBIO |Q, which is based on the WET interaction |3^. This interaction was obtained by 
fitting binding energies and energy levels, including cross-shell data. The interaction describes well the low excitation 
energy spectra for A= 10-22 nuclei. It has also been used to investigate Gamow- Teller /3 decay rates for A <18 nuclei 

H- 

In the shell model, the reduced matrix elements of transition operators are expressed as linear combinations of the 
reduced matrix elements of single-particle states with coefficients given by one-body density matrix elements 

{jfTf\\d'%m - (22) 

The form factors in (^) have to be corrected for the center-of-mass motion. This is done by multiplying the matrix 
elements by the Tassie-Barker function, exp{b^q'^ /2A), with b being the oscillator length parameter ||3q|. 



IV. DISCUSSION OF THE RESULTS 



A. Theory versus experiment 

Using the above formalism we calculate results for the flux-averaged cross sections for neutrinos coming from the 
DIF of 7r+, {Vfj^, fi~)DIF, and for neutrinos coming from the DAR of /z"*", {Vf,,e~)DAR. The flux-averaged cross 
sections are obtained by taking neutrino fluxes from As far as RPA and QRPA are concerned, we discuss in 

the following results obtained by using the force SIII. We have checked that the interaction SGII gives very similar 
results. 

Let us first discuss our results for the inclusive cross sections (table |). In the Shell Model (SM) case, cross sections 
obtained both in the Os-Op-lsOd ((0-|-l-|-2)?icj) model space and in the Os-Op-lsOd-lpOf ((0-|-1-|-2-|-3)?iijj) model space are 
shown. Results in the smaller space are only given for comparison with | [T^ . Two different types of radial wavefunctions 
have been used, the Harmonic Oscillator with h = 1.64 fm (HO wf) and Hartree-Fock wavefunctions (HF wf). HO 
wavefunctions are used even though their radial behaviour is known not to give good results of ^^C{p,n) reactions 
cross sections. The reason for using them is twofold : i) the spurious center-of-mass motion is exactly substracted 
only in this case; ii) they are employed in order to show the sensitivity of the calculated cross sections to the choice 
of the radial wavefunctions. Cross sections can be indeed quite sensitive to the radial wavefunctions since going from 
HO to HF wavefunctions in the [Uf^, ^~)DIF case produces changes of about 30%, whereas in the (i/g, e_VDv4i? case 
we have variations of a few percent only (see table |). For comparison we also show the results of ref. ||l^ obtained 
by extrapolation in the Os-Op-lsOd-lpOf space with Woods-Saxon wave functions. The results obtained within the 
Skyrme RPA and QRPA are larger than in SM. The QRPA values are sligthly larger than the RPA ones for the 
inclusive cross section. This is due to the fact that in QRPA the particle-particle residual interaction (on the average 
attractive) competes with the (on the average repulsive) particle-hole residual interaction. As a consequence, the 
QRPA strength distribution is slightly shifted towards lower energies. This fact produces a few percent increase of the 
flux-averaged cross sections due to the dependence of the differential cross sections (^ on the energy and momentum 
of the outgoing lepton and also due to the flux averaging (0,^. In ref. Ilj,^, on the contrary the introduction of 
fractional occupancies shifts the strength distribution to higher energies and this gives a few percent decrease of the 
flux averaged cross section with respect to their RPA results. From table |, we see that the (i/^t, /i~)_D/i^ cross sections 
are compatible with those of ref. [|l|, but there is a ratio of about a factor 3 for the (z/g, e~)DAR case. Most of the 
disagreement is due to the bad description of the ground state to ground state transition (its contribution represents 
2/3 of the total cross section), some of it arises because of the ground state energy of ^^iV which is not close to 
the experimental value. In plf the single particle energies were fitted to the experimental values whereas our single 
particle energies are kept as they come out from our Skyrme HF or HF-BCS calculation. The single particle energies 
are particularly important to get the energy of the ^^N ground state with respect to ^^C ground state close to its 



6 



experimental value of 17.338 MeV. In fact, within RPA this transition corresponds to the IP3/2 ^ lpi/2 transition, 
and with the Skyrme forces employed in this work the energy of this transition is 13.4 MeV. 

In table || we see the importance of including configuration mixing as well as of having the energy of this transition 
close to its experimental value. This can also be seen from table [II, where we give the inclusive {h'e,e^)DAR cross 
sections calculated without the contribution of the ground state to ground state transition. In this case, our RPA 
cross section is very close to the experimental value as found in [|l2|,^ . 

We discuss now the theoretical predictions versus the experimental findings. For the {i^fj,, ^^)DIF case, we see from 
table ^ that our results within SM (with HF wf ) are close to the experimental value when the error bars are taken into 
account, sligthly overestimating it. The WS results in agree even better. On the other hand, the RPA, the QRPA 
and Continuum RPA (CRPA) with fractional occupancies overestimate the experimental value by about 50%. 
Concerning the {Vf,,e~')DAR cross sections, we see that the SM results {Os-Op-lsOd-lpOf with HF wf) are compatible 
with the experimental values when error bars are taken into account, as those of ref. ||l^. Again the CRPA results 
of ref. overestimate the experimental values by about 30% (while our RPA predictions are far off for the reasons 
explained above). 

Let us discuss the exclusive cross sections together with various related processes like /? decay and muon capture. 
First, we can see from table || that our shell model results nicely agree with the experimental values for both the 
DIF and the DAR cases. The calculated values are not very sensitive to the choice of the radial wave functions 
and confirm those calculated in pa . Our RPA prediction for the {v^, ^~)DIF is about a factor 3 larger than the 
experimental value as found also injllj] , whereas the introduction of fractional occupancies brings the RPA predictions 
close to the experimental value [ p^21| ]. As far as the QRPA is concerned, the value obtained is slightly lower than in 
RPA, but within QRPA there are difficulties in choosing the ground state of ^"^N because the lowest state is not the 
most collective one. Therefore all the results concerning the exclusive transition have been obtained by summing the 
strength of the energy levels within the first 3 M eV above the lowest state and attributing it to a single state at an 
average energy which is not much different from the RPA one. Our RPA (z/g, e~)DAR exclusive cross section results 
are a factor 4-5 larger than the experimental values as found in A related process is the ^^iV — C (i^ decay. 
In this case the transition operator does not have a radial dependence. From table |l^, we see that the SM ft value is 
in reasonable agreement with the experimental one whereas RPA gives values four times smaller as it is well known 
0,|ll. The ft value for the /?" decay from ^"^B to is not given, but similar arguments hold since ^^B and ^"^N are 
mirror nuclei. Finally, let us look at the exclusive muon capture rates (table ^). The SM results are very close to the 
experimental value in agreement with what was found in p6| while as in pl[| RPA overestimates it by about a factor 
4. The nice agreement between experiment and theory obtained within SM (at variance with standard RPA) shows 
clearly that the inclusion of configuration mixing is necessary for a good description of the mentioned processes. 

We now turn to the discussion of the inclusive {ve, e^)DAR cross section without the contribution of the ground 
state to ground state transition (see table III). There are two main differences with the results in |l^ : i) our 
{I'e, e~)DAR cross sections are not very sensitive to the choice of the wavef unctions; ii) the calculated value is 
twice the one obtained in Jl6t . In order to understand point i), we show the strength distribution obtained for the 
multipolarity J"^ = 1+, calculated either with HO or HF wavefunctions, with q = 0.2 fm~^ in (||) as a "typical" 
value of the momentum transferred in the {i'e,e~)DAR reaction (fig.l) and q = 1.0 fm~^ as a "typical" value for 
the {i/fi, fi^)DIF case (fig. 2). We see that the strength distributions obtained for q = 0.2 fmr^ are the same for the 
two s.p. wavefunctions, while for q = 1.0 fm~^ there are significant differences. This explains the sensitivity to the 
choice of the wavefunctions of the results of table | and IH . Concerning point ii) , from (^) we see that the differential 
cross sections scale as the square of the lepton energy in the case of electrons (because of its negligeable mass) . As a 
consequence, small shifts in the strength distributions can have large effects on the cross sections as we have already 
seen in table |l|. Because neutrino-nucleus reactions with electronic neutrinos are so sensitive to the details of the 
strength distribution, they could eventually be used as a probe of nuclear structure. 

Inclusive muon capture rates, calculated according to (|9|)-(^2|) are given in table VI. The SM results in the large 
space and with HF wavefunctions are about 20% lower than in ||l^ and 10% lower than the experimental values. 
Concerning the results obtained within RPA, the disagreement between calculations (present results and ref. pi|) and 
the experimental value is again due to the bad description of the ground states, as it can be seen from table^ 

The contributions of the most important mul tipol ariti es wi th J < 6 included in the inclusive (ly^^, fi^)DIF and 
{iye,e^)F)AR cross sections are shown in tables VII and VIH respectively. We see that the contributions in RPA 
and QRPA are essentially the same, showing that the configuration mixing induced by the pairing correlations is 
not affecting very much the inclusive cross sections. In table VII the results of the SM calculations again depend 
significantly on the choice of the s.p. wavefunctions. Significant differences are also evident going from RPA to SM 
(HF wavefunctions), for example for the l"*", 1^, 2+ and 3~. On the contrary, the contributions of the different 
multipolarities to the (i'e, e~)DAR cross sections are almost equal in the three approaches, except for the remarkable 
difference associated with the 1+ states (and in particular with the ground state to ground state contribution). 

In view of future experiments using ^^C as detector for neutrinos, with impinging neutrino fluxes different from the 
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ones used up to now, in figs. 3 and 4 we give the differential cross sections obtained botli in RPA and in SM for tfie 
reactions ^^C(z^e, e~)^^iV and ^^C(;^^, /i~)^^A^ as a function of neutrino energies up to 300 MeV . 

A question arises concerning tlie reason for the difference between the RPA and SM resuhs for the inclusive 
cross section in the {v^, ijl~)DIF case (tables J,V1]), the value of which has important implications on the recent 
measurements of neutrino oscillations. The two approaches have two major differences : i) the type of correlations 
included; ii) the model space used. Concerning point i), we have already seen in the previous section th at co rrelations 
are actually responsible for the quenching in the 1+ (ground state to ground state) contribution (table VII) which is 
taken into account more correctly in the SM. Due to the tensor interaction the effects of quenching and fragmentation 
of the strength are also important for the spin-dipole mode, especially for the 1~ component (see the discussion in 
ref. [^). The difference of cross sections between SM and RPA in the case of the 0~ and 2~ is rather small, while 



it is about 0.5 x 10-*° 



for the 1 . The difference in the type of the correlations included explains half of the 



discrepancy between the SM and RPA results. Concerning ii), one important constraint on nuclear models is due to 
sum rules |39j . The existence of sum rules imposes constraints on the model spaces to be used in a given approach 
because one should ensure that the model space used is large enough to satisfy the corresponding sum rule for a given 
operator. This has not been considered in previous studies 16 

the RPA and SM cross sections. In the next subsection we discuss the role of sum rules. 



2l[| and can explain some of the difference between 



B. Sum rules 

The processes studied involve operators of two kinds, either of Fermi type ^ or of Gamow- Teller type (||). In the 
low-g limit these operators become the standard multipole operators ( p^ . As an example we have calculated sum rules 
for the non spin-flip states J'^ = 1^ and J'^ ~ 2+ both in RPA and in SM. We have found that if for the 1~ states, the 
corresponding sum rule is satisfied in the two approaches, the sum rule for the 2+ states is satisfied in the RPA whereas 
20% is missing in SM when the largest space [Os-Op-lsOd-lpOf) is used. Figure 5 shows a comparison between the 
RPA and SM strength distributions obtained for the 2+ states with the Fermi type operator, O = rf,Y2{rk)t+{k)- 
We can see that in the SM calculation some strength is missing at about 40 — 60 MeV because of the truncation of 
the space. This missing strength at high energy gives a significant contribution to the flux averaged {vf^, ^^)DIF 
cross section. 

In order to show this we have performed an RPA calculation of this cross section in the same space as used in 
the SM calculations {Os-Op-lsOd-lpOfj. The results obtained in RPA, RPA in the restricted {3Tluj) model space and 
SM are shown in table We see that for the 1~ states, for which the basis is large enough for the sum rule to 
be satisfied, the corresponding cross sections are practically the same in the three cases. On the contrary, the 2+ 
contributions to the total cross section look quite different, showing a reduction of about 10% going from RPA to 
RPA in the restricted space and to almost 30% in the SM. This reduction is due to the missing strength as illustrated 
in figure 5 and is accompanied by sum rules that are not exhausted. 

Finally, from table |^ we also see that restricting the RPA space brings the total flux averaged cross section very 
close to the SM one. This important result seems to indicate strongly that getting theoretical predictions in the 
shell- model framework which come close to the experimental values in the inclusive {up,, iJ,~)DIF can be an artifact 
because the model spaces used for the calculations are not large enough to satisfy sum rules. When a more extended 
space is used, it is quite possible that the shell model cross section will increase approaching the RPA value and 
exceeding the experimental LSND result. However, we should also note that the increase of cross sections with 
increasing configuration space is not the same for SM and RPA. The rate of increase of cross sections is smaller for 
the SM case. This can be already seen from the difference of the rate of increase when going from ITiuj to 3huj spaces 
when we compare our SM results with ref. ||l^. The rate of increase is 14% in our case while it is 19% in ref. [ p^ , 
where the cross section for Shuj space was obtained from an extrapolation based on an RPA-like model. 



C. The problem of quenching of spin strength 

Finally, we would like to discuss more the problem of quenching of spin strength in nuclei. It is known for sometime 
that there is not enough strength in the main GT peak. Two explanations on the origin of this quenching have been 
suggested : i) the missing strength is shifted to very high energy due to the coupling of the nucleon internal structure 
to the Delta resonance; ii) there is a shift of the strength to energies up to about, or more than 50 MeV due to the 
coupling to multi-particle multi-hole configurations. The quenching and fragmentation of the spin strength due to 
these effects would lead to quenching of the cross sections of the GT mode and probably for the spin-dipole modes 
as well. We have made a test calculation in which we have introduced an effective axial vector coupling constant. 
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(7^ . Considering high energies of neutrinos in the DIF experiments most of the missing strength due to the above 
multiparticle effects is picked up and also some portion of the strength due to the couphng to the Delta-Resonance 

. Therefore the g^l^ one should use here should be closer to qa as compared to the value 
found for the main GT peak {jfj^^ = 0.85^). 

In table we show all the results concerning the different processes using a quenching factor. For the sake of 
demonstration we use g"]^^ — 0.9gA- Using this value, the /i value is increased as expected by 20% and the exclusive 
and inclusive muon capture is reduced by about 15 — 20%. Concerning the exclusive cross sections are reduced by 
about 15% and become smaller than the observed ones, whereas the inclusive ones are reduced by about 10 — 15% 
and get closer to the experimental values. The exclusive neutrino cross sections we have seen are very sensitive to the 
calculated wavefunctions and therefore cannot be a good test of the quenching factor. The inclusive cross sections 
are less sensitive to nuclear structure and they seem to favor some reduction in qa- There are many uncertainties in 
these arguments. We do not yet know for sure what is the mechanism of quenching of qa (although the experiment 
in [ ^sf suggests that the main role is played by the coupling with many particle-many hole states) and whether the 
same quenching factor should be applied to all multipolarities. It is quite conceivable that future studies of the 
neutrino-nucleus interaction will actually help to clarify this point about quenching. 



effect is recovered |40 



V. CONCLUSIONS 



Microscopic approaches, namely charge-exhange RPA, charge-exhange QRPA and Shell Model, are used to evaluate 
^'^C{yf,,e~Y^N and ^^C(i'^, [T'Y'^N cross sections both for Vf. coming from the decay-at-rest of /x"*" and for coming 
from the decay- in- flight of 7r+ . Accurate knowledge of these cross sections is important for the interpretation of recent 
measurements of neutrino oscillations performed both by LSND and KARMEN collaborations and also for future 
experiments. The results show that the calculated exclusive cross sections, where ^"^N is left in the ground state, are 
in good agreement with the measured values when large-scale shell model calculations are performed. In fact, in this 
framework it is easy to include properly the configuration mixing present in the ground state of ^^C and therefore to 
have a good description of the ground state wavefunctions. Concerning the inclusive ^^C(i/^, /i~)^^7V cross sections 
with from DIF of 7r+, in which ^^A^ is left either in the ground state or in an excited state, we get in the shell model 
15.2 X 10^'*" err? (when Hartree-Fock wavefunctions are used), about 20% larger than the experimental value, which 
is 12.4 ± 0.3 ± 1.8 X 10"''*' err? . The calculated cross section in charge-exchange RPA is 50% larger than this value, 
namely 19.2 x 10~'*° err? . The most important difference between the SM and RPA is the quenching in the 1+ due to 
the different correlations included in the two approaches. This explains half of the difference between the two results. 
The evaluation of sum rules for natural parity states has also shown that the basis used in the shell model calculations 
is not large enough to exhaust these sum rules whereas in the case of the RPA the sum rules are satisfied. This fact 
suggests that the reduced cross section obtained in the shell model framework is partially due to the use of a basis 
which is not large enough. Enlarging the model space would then add some strength at high energy and therefore 
increase the inclusive cross sections. From these arguments, we conclude that in both microscopic approaches the 
theoretical prediction is 20 — 30% larger than the measured value. But even if we could hope to extend further the 
shell model calculation and increase the basis, the calculated value will keep having 10 — 20% uncertainty due to a 
certain degree of arbitrariness in the choice of the wavefunctions for the unbound states and also of the interactions. 
It would be interesting to extend the shell model calculations to larger configuration space in spite of the uncertainties 
above, to have more information on the correlations in spin-dipole modes and clarify the convergence of the cross 
sections as the space is extended. Concerning the inclusive ^'^Ciy^.e^Y'^ N cross section we get 16.4 x 10~^^ err? 
which agrees within the experimental error bars with the measured value. The ft value for the beta decay from ^^iV 
to ^^C, exclusive and inclusive muon capture rates ^"^Bi^pT .v^^'^C are also evaluated and are in quite good agreement 
with the measured values. 

This work was supported by the US-Israel Binational Science Foundation. The shell model calculation was done 
by using the TKYNT at Department of Physics, University of Tokyo. 
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TABLE I. Flux averaged inclusive cross sections < a >f within the different approaches used, namely the Shell Model (SM), 
the Random Phase Approximation (RPA) and the Quasi-particle RPA (QRPA). Different SM results are given according to 
various choices of radial wavefunctions, i.e., oscillator functions with length parameter b = 1.64 fm (HO wf) and Hartree-Fock 
wave functions (HF wf) and for either (0 + 1 + 2)hu! or (0 + 1 + 2 + 3)hLj model space. Comparison with CRPA with fractional 
occupancies and a recent shell model calculation ||l^ is made. In the former case the results are obtained with the 
finite-range G- matrix derived from the Bonn NN potential (BP) and with the Landau-Migdal (LM) force (in brackets). In the 
latter case, Woods-Saxon wavefunctions (WS wf ) have been used and the results within the (0 -|- 1 -|- 2 -|- 3)huj model space (in 
brackets) are obtained by extrapolation. 



(i.^,At-)D/F 


(i^e,e-)DAR 


< a >f (10-«' cm^) 


<(T>f {IQ-^^ cm^) 



SM(HO wf) (0 -f 1 + 2)huj 

SM(HF wf) (0 + 1 + 2)hLo 

SM(WS wf) (0 -f 1 -I- 2)hLo ^ 

SM(HO wf) (0 + 1 + 2 + ?,)huj 

SM(HF wf) (0 + 1 + 2 + ?,)hLj 

SM(WS wf) (0 + 1 + 2 + ?.)hLo 

RPA 

QRPA 

CRPA ll) 

EXP 







18.71 
13.33 
11.1 
21.08 
15.18 
(13.2) 
19.23 
20.29 
18.18(17.80^ 
12.4 ± 0.3 ± 1.8 



§1 



14.21 
13.94 
12.1 
16.70 
16.42 
(12.3) 
55.10 
52.0 
19.28(18.15) 
14.1 ± 1.6 ± 1.9 
14.8 ± 0.7 ± 1.4 
14.0 ±1.2 § 



TABLE II. Same as table 1 for flux averaged exclusive cross sections. 





{u,,,ti~)DIF 


{iye,e')DAR 




<a>f (10-*" cm^) 


<a>f {10-^^ cm') 


SM(HO wf) + 2)huj 


0.70 


8.42 


SM(HF wf) {0 + 1 + 2)huj 


0.65 


8.11 


SM(WS wf) (0 + 1 + 2)hui ^ 


0.58 


8.4 


RPA 


2.09 


49.47 


QRPA 


1.97 


42.92 


CRPA (21) 


1.06(1.03) 


13.88(12.55) 


EXP 


0.66 ± 1.0 ± 1.0 § 


10.5 ± 1.0 ± 1.0 jfl 






9.1 ±0.4 ±0.9 m 






9.1 ±0.5 ±0.8 gl 
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TABLE III. Same as table 1 for flux averaged inclusive cross sections but excluding the ground state. 



{u^,e~)DAR 

SM(HO wf) (0 + 1 + 2)hLo 5.79 
SM(HF wf) (0 + 1 + 2)huj 5.83 
SM(WS wf) (0 + 1 + 2)nuj H 3.7 
SM(HF wf) (0 + 1 + 2 + ■i)hu 8.28 
SM(HO wf) (0 + 1 + 2 + 'i)huj 8.31 
SM(WS wf) (0 + 1 + 2 + 3)^1.^ 10 (3.8) 
RPA 5.63 
QRPA 9.08 
CRPA (21) 5.4(5.6) 
EXP 5.4 ± 1.9 1] 

5.7 ±0.6 ±0.6 § 
5.1 ±0.8 § 



TABLE IV. ft value for the /3+-decay from ^^Ngs to ^^Cgs 



SM(HF wf) (0 + 1 + 2)huj 17008 
SM(HO wf) (0 ± 1 + 2)hio 16425 
RPA 3032 

EXP 13182 ± 1. @ 



TABLE V. Same as table 1 for the exclusive muon capture rates Ac. 



(10-* s~^) 

SM(HO wf) (0 + 1 + 2)hLj 0.50 

SM(HF wf) (0 + 1 + 2)huj 0.48 

SM(WS wf) (0 + 1 + 2)hij ^ 0.66 

RPA 2.54 

CRPA (21) 2.37(2.43) 

EXP 0.62 ± 0.03 H 



TABLE VI. Same as table 1 for the inclusive muon capture rates Ac. 



(10-* s-') 

SM(HF wf) (0 + 1 + 2 + 3)hij 3.32 
SM(WS wf) (0 + 1 + 2 + 3)a(x) ||l§ (4.06) 
RPA 5.12 
CRPA |2l) 5.79(5.76) 
EXP 3.8 ±0.1 @ 
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TABLE VII. Contribution of the most important multipolarities to the inclusive {v^, n~)DIF cross sections 
< a >f (10"*° cm^), within the RPA, QRPA and SM approaches. For the latter, results obtained with two different choices 
of the wavefunctions, Hartree-Fock (HF wf) and Harmonic Oscillator (HO wf) are shown. 



J" 


RPA 


QRPA 


SM 


SM 








HF wf 


HO wf 


0" 


0.96 


0.95 


0.82 


1.35 


1+ 


4.42 


4.48 


2.47 


4.11 


1" 


3.53 


3.48 


3.11 


4.06 


2+ 


2.04 


1.94 


1.38 


2.29 


2" 


3.78 


4.21 


3.87 


4.53 


3+ 


1.71 


2.25 


1.58 


2.33 


3" 


0.70 


0.80 


0.47 


0.58 


4" 


1.36 


1.35 


1.11 


1.19 



TABLE VIII. Contribution of the most important multipolarities J" to the inclusive {vs,e~)DAR cross sections 
< a >f (10~*^ cm^), within the RPA, QRPA and SM approaches. For the latter, results obtained with two different choices 
of the wavefunctions, Hartree-Fock (HF wf) and Harmonic Oscillator (HO wf) are shown. 





RPA 


QRPA 


SM 


SM 








HF wf 


HO wf 


0" 


0.3 


0.5 


0.6 


0.6 


1+ 


50.0 


45.6 


9.6 


9.8 


1" 


1.7 


2.0 


2.1 


2.2 


2+ 


0.1 


0.1 


0.1 


0.1 


2" 


3.0 


3.7 


4.0 


3.9 


3+ 


0.0 


0.1 


0.0 


0.1 


TABLE IX. Comparison of the the inclusive {Vfi,, 


')DIF cross sections < cr >/ (10 *° cm^) obtained within RPA, RPA in 


the same model space as the SM and SM with HF wavefunctions. 






J'^ 


RPA 




RPA(3/iu;) 


SM (HF wf) 


0+ 


0.15 




0.14 


0.15 


0" 


0.96 




0.81 


0.82 


1+ 


4.42 




3.76 


2.47 


1" 


3.53 




3.14 


3.11 


2+ 


2.04 




1.78 


1.38 


2" 


3.78 




3.09 


3.87 


3+ 


1.71 




1.23 


1.58 


3" 


0.70 




0.40 


0.47 


4+ 


0.18 




0.07 


0.07 


4" 


1.36 




1.13 


1.11 


5+ 


0.39 




0.24 


0.15 


6+ 


0.01 




0.00 


0.00 


TOTAL 


19.23 




15.79 


15.18 
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TABLE X. Results for the different processes obtained with g'J^' = 0.9gA both within RPA and SM. 



SM (HF wf) 

{v^,H-)DIF <a>f (10-™ crri^) exclusive L76 055 

inclusive 17.02 13.49 

{ve,e-)DAR<a>f {IQ-'^'^ cm^) exclusive 40.79 6.69 

inclusive 45.63 13.8 

ft value(s) 3743 20278 

Ac (10'' s"^) exclusive 2.05 0.387 

inclusive 4.11 2.78 
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FIG. 1. Strength distributions for the J'^ = 1"*" states obtained with the transition operator 6 = t+{k)j(,{qrk)\Y(,{rk) xa] 
with g = 0.2 fm~^, using both HO (full line) and HF wavefunctions (dotted line). The lines are the results of a folding with 
lorentzian of 1 MeV width. 
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FIG. 2. Strength distributions for the J'^ — 1'^ states obtained with the transition operator O — t+{k)j(){qrk)[Y(){rk) x a] 
with q = 1.0 fm~^ , using both HO (full line) and HF wavefunctions (dotted line). The lines are the results of a folding with 
lorentzian of 1 MeV width. 



15 




E (MeV) 

FIG. 3. Inclusive {ue,e^)DIF differential cross section as a function of neutrino energy, calculated both in RPA (full line) 
and in SM (dotted line). 




E (MeV) 

FIG. 4. Inclusive {v^, IJ,~)DIF differential cross section as a function of neutrino energy, calculated both in RPA (full line) 
and in SM (dotted line). 
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FIG. 5. Strength distributions for the = 2+ obtained with the Fermi type operator O = r^Y2{rk)t+{k) with the RPA 
approach (full line) and SM (dotted line) . The lines are the results of a folding with a lorentzian of 1 MeV width. 
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